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Abstract 

Using the low-energy effective field theory for magnons and holes 
— the condensed matter analog of baryon chiral perturbation theory for 
pions and nucleons in QCD — we study different phases of doped antifer- 
romagnets. We systematically investigate configurations of the staggered 
magnetization that provide a constant background field for doped holes. 
The most general configuration of this type is either constant itself or it 
represents a spiral in the staggered magnetization. Depending on the val- 
ues of the low-energy parameters, a homogeneous phase, a spiral phase, 
or an inhomogeneous phase is energetically favored. The reduction of 
the staggered magnetization upon doping is also investigated. 
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1 Introduction 



The precursors of high-temperature superconductors [Jj are doped antiferromagnets 
with a spontaneously broken global SU (2)^ spin symmetry and with magnons as the 
corresponding Goldstone bosons. The effect of antiferromagnetic spin fluctuations 
on the dynamics of doped holes has been investigated in great detail in the condensed 
matter hterature BJpJi BlBB iiBEHJll^ ■ 

H H H H |23,M|2i liila,!!!^^ 

Using a variety of numerical and analytic techniques, a wide range of interesting 
phenomena has been investigated in doped antiferromagnets. In particular, it was 
suggested that spiral phases with an inhomogeneous staggered magnetization may 
replace the Neel pha se of the undoped antiferromagnet even at arbitrarily small 

(1 Q El ll 111 H 111 lil H H 111, H H H H H in 



doping 

a spiral phase the staggered magnetization develops a helix structure, and the Neel 
ordered antiferromagnet thus turns into a helimagnet. Other inhomogeneities — 
most important stripes — have also attracted a lot of attention ji^ . Unfortunately, 
away from half-filling, the microscopic Hubbard and t-J models cannot be simulated 
reliably due to a severe fermion sign problem. Also analytic calculations are usually 
not fully systematic but suffer from uncontrolled approximations. Consequently, 
most results for these strongly correlated systems remain, at least to some extent, 
debatable. While this may seem unavoidable taking into account the complicated 
nonperturbative dynamics of these systems, a systematic effective field theory ap- 
proach allows us to reach some unambiguous conclusions at least for lightly doped 
systems. While some results of this paper have been derived before using less rig- 
orous methods, the effective field theory derivation is still very useful, because it is 
reliable and particularly transparent. 

Particle physicists are facing the challenges of strongly correlated systems in 
studies of the strong interactions between quarks and gluons. Just like an undoped 
antiferromagnet, the QCD vacuum has a spontaneously broken global symmetry 
— in that case the SU{2)l ® SU{2)ii chiral symmetry — which gives rise to three 
Goldstone pions — the analogs of the magnons in an antiferromagnet. The QCD 
analog of the doped holes carrying electric charge are the nucleons carrying baryon 
number. Just as simulating the Hubbard model at non-zero doping is prevented by 
a fermion sign problem, simulating QCD at non-zero baryon density is prevented by 
a severe complex action problem. For this reason, lattice QCD is presently limited 
to simulating individual particles propagating in the QCD vacuum. Although sim- 
ulations of the QCD vacuum do not suffer from the complex action problem, they 
are still very demanding, especially in the physical regime of small quark masses. 
Fortunately, a systematic effective field theory ^ il 47, 4^ — chiral perturbation 
theory — is extremely successful in describing the low-energy physics in this regime. 
In chiral perturbation theory not quarks and gluons but pions and nucleons are the 
fundamental degrees of freedom. Although the effective theory is not renormalizable, 
it allows a systematic low-energy expansion with only a finite number of a priori un- 
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known low-energy parameters entering at each order. The values of the low-energy 
parameters can be determined from experiments or from lattice QCD simulations. 
Chiral perturbation theory provides us with precise predictions for the low-energy 
pion physics, which would be practically irapqssible to derive directly from QCD. 
Baryon chiral perturbation theory ji^ 0, IHH |H3, extends these successes to 
the low-energy physics of both pions and nucleons. At present, a fully systematic 
power-counting scheme seems to exist only for the sector with a single nucleon 
Still, few-nucleon systems have also been treated quantitatively j5J, 



531. 



The QCD analog of a smral phase in a doped antiferromagnet is a pion condensate 
in nuclear matter 0, S H E H • 

The systematic technique of chiral perturbation theory is not limited to QCD 
but can be applied to any system with Goldstone bosons. Indeed, systematic low- 
energy effective theories have been very successful in describing the dynamics of 



ma gno ns in both ferro- and antiferromagnets j9|, llO|, ll3|, |23|, |30|, |35|, |65|, |66|, |67, 

we have extended the pure magnon effective field theory by 
including charge carriers. The resulting effective theory for magnons and doped 
holes is the condensed matter analog of baryon chiral perturbation theory. The 
effective theory incorporates important experimental as well as theoretical results, 
such as the location of hole pockets at lattice momenta (^, which follows from 
angle resolved photoemission spectroscopy (ARPES) experiments [22, .7 3.. 7 4 as 
well as from theoretical investigations of Hubbard or t-J-like models ja, [D, 120, |4l|. 
Recently, we have used the effective theory to derive the m agno n-mediated forces 
between two isolated holes in an otherwise undoped system [tiI . Fg^ . Remarkably, 
the Schrodinger equation corresponding to the one-magnon exchange potential can 
be solved analytically and gives rise to an infinite number of bound states. It remains 
to be seen if these isolated hole pairs are related to the preformed Cooper pairs of 
high-temperature superconductivity. In this paper, we use the effective theory to 
investigate the regime of small doping. This is possible analytically if the 4-fermion 
contact interactions between doped holes are weak and can be treated perturbatively. 
Whether this is indeed the case depends on the concrete magnetic material under 
consideration. It should be noted that the 4-fermion couplings between doped holes 
in the effective theory may well be small, although the microscopic on-site Coulomb 
repulsion U in the Hubbard model or the exchange coupling J in the t-J model 
which cause antiferromagnet ism are strong. In particular, in the effective theory 
antiferromagnetism arises independent of the strength of the remnant 4-fermion 
couplings between doped holes. Assuming that the 4-fermion couplings can be 
treated perturbatively, the effective theory predicts both homogeneous and spiral 
phases, depending on the specific values of the low-energy parameters. 

The paper is organized as follows. In section 2 the effective theory of magnons 
and holes as well as the nonlinear realization of the spontaneously broken SU{2)s 
spin symmetry are reviewed. The holes interact with the Goldstone bosons via a 
U{l)s "gauge" field and two "charged" vector fields composed of magnons. The 
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gauge group U{l)s and the corresponding "charge" refer to the unbroken subgroup 
of SU{2)s. In section 3 we consider the most general magnon field that gives rise to 
constant gauge and charged vector fields and thus to a homogeneous background for 
the doped holes. These magnon fields turn out either to be homogeneous themselves 
or to form a spiral in the staggered magnetization. A particular magnon field which 
gives rise to inhomogeneous composite gauge and charged vector fields — a so- 
called double spiral — is also discussed. In section 4 homogeneous and spiral phases 
are investigated. The effect of weak 4-fermion contact interactions is investigated in 
section 5 using perturbation theory and — depending on the values of the low-energy 
parameters — it is determined which phase is energetically favored. In section 6 
the reduction of the staggered magnetization upon doping is calculated both for 
the homogeneous and spiral phases. Section 7 contains an outlook as well as our 
conclusions. In an appendix we prove that the most general configuration of the 
staggered magnetization that provides a constant background field for the doped 
holes is either constant itself or represents a spiral. 



2 Systematic Low-Energy Effective Field Theory 
for Magnons and Holes 



In order to make this paper self-contained, in this section we review the effective 
theory for magnons and holes constructed in 7(1, 7]]| which is based on the pure 
magnon effective theory of |l3 Q HI, H IH, IfillilB Q . 



The staggered magnetization of an antiferromagnet is described by a unit-vector 



field 



eix) 



(sin6'(x) cos(p{x), sin6'(a;) sm(p{x), cos9{x)), (2.1) 

in the coset space SU {2)s/U{l)s = S^, with x = (xi, X2, t) denoting a point in (2+1)- 
dimensional space-time. For our purposes it is more convenient (but completely 
equivalent) to use a CP(1) representation in terms of 2 x 2 Hermitean projection 
matrices P{x) that obey 



P{x) = -{l + e{x) - a), P{x)^ = P{x), TrP(x) = 1, P{xf = P{x), (2.2) 
where a are the Pauli matrices. The relevant symmetries are realized as follows 



^^7(2), 


P{x)' 


= gP{x)g\ 










^'P(x) 


= 1 - P(x), 











°P(x) 


= P{Ox), 
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{-X2,Xi 
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^P(x) 


= P{Rx), 
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(Xi, —X2 
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^P(x) 


= l-P{Tx), 


Tx = 


{xi,X2, - 


-t). 



(2.3) 
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Here g G SU{2)s is a matrix that implements the global spin symmetry which is 
spontaneously broken down to U{l)s, Di denotes the displacement by one lattice 
spacing in the i-direction, and O, R, and T denote 90 degrees spatial rotations, 
spatial reflections, and time-reversal, respectively. 

In order to couple doped holes to the magnons, a nonlinear realization of the 
SU{2)s symmetry was constructed in [l^. The magnon field is diagonalized by a 
unitary transformation u{x) G SU{2)s, i.e. 

M(x)P(x)M(x)t = i(l + a3) = (^ J Mn(x)>0, (2.4) 

with 

u(x)=( cos{\e{x)) sin(|^(x))exp(-i(^(a;)) 

Y — sin(i^(a;)) exp(i(y9(a;)) co?,{\9{x)) 

The transformation u{x) describes a rotation of the local staggered magnetization 
vector e(x) into the 3-direction. Since u{x) is more directly related to P{x) than 
to e{x) itself, we have chosen the CP(1) representation. Under a global SU{2)s 
transformation g the diagonalizing field u{x) transforms as 

u{xy = h{x)u{x)g\ Mn(x)'>0, (2.6) 

which defines the nonlinear U{l)s symmetry transformation 

h{x) = exp{ia{x)(T3), (2.7) 

that acts like a gauge transformation in the unbroken subgroup U{l)s. The local 
symmetry transformation h{x) depends on the global transformation g as well as 
on the local staggered magnetization P{x) from which it inherits its x-dependence. 
Under the displacement symmetry Di one obtains 

""W-w). - ( -^p'-^w). (2.8) 




The way in which the global SU{2)s spin symmetry disguises itself as a local 
symmetry in the unbroken U{l)s subgroup is characteristic for any systematic ef- 
fective field theory of Goldstone bosons. The nonlinear realization of spontaneously 
broken continuous global symmetries has been discussed in full generality in the 
pioneering work of Callan, Coleman, Wess, and Zumino [i^ |4^]. Following their 
general scheme, doped holes are derivatively coupled to the magnons. In fact, the 
holes are "charged" under the local U{l)s symmetry and transform with the non- 
linear transformation h{x). In order to couple holes to the magnons it is necessary 
to introduce the anti-Hermitean traceless field 

Vi,{x) = u{x)d^u{x)\ (2.9) 
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Figure 1: Elliptically shaped hole pockets centered at (±^,±^). Two pockets cen- 
tered at and ^'^ + ? ^) combine to form the pockets for the flavors f = a, p. 



which obeys the following transformation rules 



SU{2), 




= h{x)[v^{x) 


+ d^]hixy 








= t{x)[v^{x) 











= EijVjiOx), 
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vt{0x), 
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= 


-V2{Rx), 
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^Vj{x) 


= ""'v.iTx), 


'^Vt{x) = 


- ^^VtiTx 



%t{x) = vt{Rx), 



Writing 



tv^[x)aa 



v^{x) 
)mposi 

which transform as 



(2.10) 
(2.11) 



the field Vfj,{x) decomposes into an Abelian "gauge" field v^{x) and two "charged" 
vector fields f^(a;) 



vlix)' 



df,a{x) 



v^ix)' = v^ix)exp{±2ta{x)), (2.12) 
under Sf/(2),. 

ARPES measurements jll, [ij 24 , [z!] as well as theoretical calculations in t- 



J-like models p, |20|, |41[ have revealed that at small doping holes occur in pockets 
centered at k°' = (^, ^) and fc^ = (^, — ^) in the Brillouin zone. The elliptically 
shaped hole pockets are illustrated in figure 1. The effective field theory is defined 
in the space-time continuum and the holes are described by Grassmann- valued fields 
ipl{x) carrying a "flavor" index f = a,f3 that characterizes the corresponding hole 
pocket. The index s = ± denotes spin parallel (+) or antiparallel (— ) to the local 
staggered magnetization. Following |70|,|7l|, under the various symmetry operations 
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the hole fields transform as 



SU(2). 




= exp(ibiQ;(x))'?/^£(x), 


UiDn 




— exp(ici;)'0£(x), 




^'ibi(x) 


= TexDfz/c'^ a) exD(^^i(D( x))ibL( x) 


o 
\j 


'ip±[x) 


— ^V±y^x}^ y±\X) — w±\yx}^ 


R 




^i^liRx), ''i^i{x)^r±{Rx), 


T 




^TeMTMTxWJiTx), 






= ± exp(±^^p(Txj)^p:^ (Tx) . 



(2.13) 

Here U{1)q is the fermion number symmetry of the holes. Interestingly, in the 
effective continuum theory the location of holes in lattice momentum space manifests 
itself as a "charge" k( under the displacement symmetry Di. 

Once the relevant low-energy degrees of freedom have been identified, and the 
transformation rules of the corresponding fields have been understood, the construc- 
tion of the effective action is uniquely determined. The low-energy effective action of 
magnons and holes is constructed as a derivative expansion. At low energies terms 
with a small number of derivatives dominate the dynamics. Since the holes are 
heavy nonrelativistic fermions, one time-derivative counts like two spatial deriva- 
tives. Here we limit ourselves to terms with at most one temporal or two spatial 
derivatives. One then constructs all terms consistent with the symmetries listed 
above. The effective action can be written as 



S[^l>i\^l^i,P] = J d'xdt 



(2.14) 



where denotes the number of fermion fields that the various terms contain. The 
leading terms in the pure magnon sector take the form 



Co = PsTr[diPdiP + -dtPdtP] 



Pi 
2 



die- die+ ^dte- dte 



(2.15) 



Here ps is the spin stiffness and c is the spinwave velocity. The leading terms with 
two fermion fields (containing at most one temporal or two spatial derivatives) are 
given by 



E 



1 



+ 



2M' '2M' 



(2.16) 
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It should be noted that vf{x) contains a spatial derivative, such that magnons and 
holes are indeed derivatively coupled. In eq. ()2.16p M is the rest mass and M' and 
M" are the kinetic masses of a hole, A is a hole-one-magnon, and A^^i and N2 are 
hole-two- magnon couplings, which all take real values. The sign o"/ is + for / = a 
and — for f = {3. The covariant derivative takes the form 



The leading terms with four fermion fields are given by 



(2.17) 



s'=+ - 
s'=+- 



(2.18) 



with the real- valued 4-fermion coupling constants Gi, G2, G3, and G4. Here we have 
limited ourselves to terms containing at most one spatial derivative. In principle 
there are even more contact interactions among the fermions such as 6- and 8- 
fermion couplings as well as 4-fermion couplings including more derivatives. Some 
of these terms have been constructed in ||2l| but play no role in the present work 
and have hence been suppressed. 

Remarkably, the above Lagrangian has an accidental global U{1)f flavor sym- 
metry that acts as 

U{1)f: ^ipiix) =exp{afiri)tlji{x). (2.19) 

This symmetry is not present in the underlying microscopic systems and is indeed 
explicitly broken by higher-order terms in the effective action. For c — > 00 the 
leading terms of the effective action also have an accidental Galilean boost symmetry 



G : 



^P{x) = P{Gx), Gx = {x-vt,t), 
^ilj:j.{x) = exp (ip-^ ■ x — ujH^ ipj^{Gx), 
(x) 



il)^^{Gx) exp i^—ip^ ■ X + 00^ t) 



(2.20) 



with p^ = {p{,pI) and cj-^ given by 
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Vl 
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(2.21) 
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Also the Galilean boost symmetry is explicitly broken at higher orders of the deriva- 
tive expansion. In the real materials Galilean (or actually Poincare) invariance is 
spontaneously broken by the formation of a crystal lattice, with phonons as the 
corresponding Goldstone bosons. Here we assume that phonons play no major role 
in the cuprates, and we focus entirely on the magnons. Still, phonons and a spon- 
taneously broken Galilean symmetry could be included in the effective field theory 
if necessary. 



3 Spirals in the Staggered Magnetization 

In the following we will consider configurations e{x) of the staggered magnetiza- 
tion which — although not necessarily constant themselves — provide a constant 
background field for the doped holes. We restrict ourselves to time-independent 
configurations, such that Vt{x) = 0. The most general configuration, with Vi{x) 
constant up to a gauge transformation, represents a spiral in the staggered magne- 
tization. We also discuss a so-called double spiral which gives rise to a non-uniform 
composite vector field and thus to an inhomogeneous fermion density. 

3.1 Spirals with Uniform Composite Vector Fields 

Since the holes couple to the composite vector field Vi{x) in a gauge covariant way, 
it is sufficient to assume that Vi{x) is constant only up to a gauge transformation, 
i.e. 

9(x) 

"^ii^y = '^ii^) " dia{x) = sin^ diip{x) — dia{x) = c^, 

v^{xy = (x) exp{±2ia{x)) 

— l-[sm9{x)diip{x) ±idi9{x)\ exp{'^i{if{x) —2a{x))) — cf, (3.1) 

with cf and cf being constant. As shown in the appendix, the most general configu- 
ration that leads to a constant Vi{xy represents a spiral in the staggered magnetiza- 
tion. Furthermore, by an appropriate gauge transformation one can always achieve 

4 = cr = c, e M. (3.2) 
The magnon contribution to the energy density of these configurations is given by 

Cm = ^di^{x) • die{x) = 2psvf{x)v~{x) = 2ps{cl + cj). (3.3) 

To be specific, let us consider a concrete family of spiral configurations with 

9{x) = 9o, (p{x) = kiXi, (3.4) 
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which imphes 

vt{x) = 0, vf{x) = kism^—, vf{x) = -^sm9oexp{TikiXi). (3.5) 
Choosing the gauge transformation 

a{x) = ^kiXi, (3.6) 

one obtains the constant field 

vtix)' = 0, vfix)' = vfix) - 9,«(x) = h{sm' y - ^) = cl 

k- 

vf{x)' = vf{x)exp(±2ia{x)) = -^smOo = Ci. (3.7) 
Hence, comparing with the appendix, in this case we can identify 

= — cos^o, a = ^ = — tan^o, (3.8) 
and the magnon contribution to the energy density is 

= 2p,(c? + 4) = y (A;? + k^) sin^ Oo- (3.9) 

3.2 A Double Spiral 

For most of this paper we restrict ourselves to configurations of the staggered mag- 
netization which give rise to a homogeneous composite vector field fi(x)'. However, 
in this subsection we examine a configuration with an inhomogeneous Vi{x)' — the 
so-called double spiral [l^. Although we will not explore this configuration any 
further in this work, it is interesting for future investigations. In particular, one can 
study it in order to figure out whether spirals with constant Vi{x)' may be unstable 
against developing inhomogeneities. 

A particularly simple form of a double spiral is given by 

e(x) = (sin(A;iXi) cos(A;2a;2), sin(A;2a;2), cos(A;iXi) cos(A;2X2)). (3.10) 

The magnon contribution to the energy density of the double spiral takes the form 

em = y<9ie(x) ■ die{x) = y [kf cofi^{k2X2) + kf) . (3.11) 
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It is straightforward to compute the composite vector field for the double spiral and 
one obtains 



VI{X) : 



ki cos(A;ia;i) sin(A;2a;2) cos(A;2.X2) 
2(1 + cos(A;iXi) cos(A;2X2)) ' 
ki cos(/c2a;2)[sin(/cia;i) sin(/c2a;2) ± i{cos{kiXi) + cos(A;2a;2))] 
2(1 + cos(A;ia;i) cos(/c2a;2)) 
k2 sin(A;iXi) 



2(1 + cos(A;iXi) cos(A;22^2)) ' 



±. _ fe2[cos(A:iXi) + cos(fc23:2) T i sin(A:iXi) sin(A:2a:2)] 
2(l + cos(A;ia;i)cos(A;22;2)) ' 

As before we perform a gauge transformation 

(x)' = vf{x) - dia{x), vfix)' = vf{x) exp{±2ia{x)) , (3.13) 
in this case with 

ex.D(2ia(x)) = ^^^(^i^i) sin(to) - ^(cos(fci3:i) + cos(fc2a:2)) , 

^ 1 + cos(A;ia;i) cos(A;2a;2) ' ^ ' ' 

and we obtain the remarkably simple form 

ki ki 
^1(2^)' = -y sin(/c2a;2), vf{x)' = — co^{k2X-2), 

^2^ = 0, vt{x)'^Ti^. (3.15) 



4 Homogeneous versus Spiral Phases 



In this section we calculate the fermionic contribution to the energy density of a 
configuration with constant Vi{x)' in order to decide which configuration is energet- 
ically favored. It will turn out that this depends on the values of the low-energy 
parameters. For large ps the magnon contribution to the energy density dominates, 
and a homogeneous phase is favored. In that case, all four hole pockets are equally 
populated with doped holes of both spin up and spin down. For smaller p^, on the 
other hand, the fermionic contribution to the energy density dominates and favors 
a spiral configuration. In this case, only a particular linear combination of spin up 
and spin down states is occupied by doped holes. If the spiral is oriented along a 
crystal lattice axis (a zero degree spiral) hole pockets of both types [a and /5) are 
filled with doped holes. On the other hand, if the spiral is oriented along a lattice 
diagonal (a 45 degrees spiral) either three or one hole pocket arc populated. As we 
will see, the zero degree spiral is realized for intermediate values of p<j, while the 45 
degrees spiral is unstable against the formation of inhomogeneities in the fermion 
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density. Interesting calculations of a similar nature were presented in j42l . |43| in 
the context of microscopic t-J-like models. In these works, in a particular parame- 
ter range a zero degree spiral has been identified as the most stable configuration. 
Our effective field theory treatment complements these results in an interesting way. 
First, it is model-independent and thus applicable to a large variety of microscopic 
systems, because the most general form of the effective action is taken into account. 
In addition, it is controlled by a systematic low-energy expansion. Before we discuss 
homogeneous versus spiral phases, we address the issue of phase separation in the 
context of the t-J model. 



4.1 Stability against Phase Separation 

It is well-known that the t-J model shows phase separation for small values of t. In 
this case the doped holes are heavy. Each hole is surrounded by four bonds which, 
in the absence of the hole, would carry a negative antiferromagnetic contribution 
to the energy. In order to minimize the number of broken antiferromagnetic bonds, 
the holes may accumulate in some region of the lattice, thus leaving an otherwise 
undoped antiferromagnet behind, i.e. the system undergoes phase separation. 

The energy density of the undoped antiferromagnet is eo. In the t-J model, which 
reduces to the Heisenberg model at half-filling, the energy density was determined 
with a very efficient loop-cluster algorithm as eo = —0.6693(1) J/a^, where J is the 
exchange coupling and a is the lattice spacing [l^. A doped hole propagating in 
the antiferromagnet has mass M. Hence, to leading order in the fermion density n, 
the energy density of a doped antiferromagnet is 

e = eo + Mn + C(n2). (4.1) 

When the system phase separates, the doped holes accumulate in some region of 
volume V , leaving an otherwise undoped antiferromagnet behind in the remaining 
volume V — V. If there are no electrons at all in the hole-rich region, each hole 
occupies an area and hence V = Na'^, where N is the number of holes. In the 
t-J model, the unoccupied region of volume V does not contribute to the energy 
density, while in the remaining volume V — V one has the energy density eo of an 
undoped antiferromagnet. Hence, the energy density of a phase separated system is 

e,s = ''^^~^'^ =eo{l-a'n). (4.2) 

The doped t-J model is thus stable against phase separation as long as 

e < eps M < -eoa^ = 0.6693(1) J. (4.3) 

In the t-J model at small t the holes are essentially static, breaking four antifer- 
romagnetic bonds, while eo represents the energy of two intact antiferromagnetic 
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bonds. Hence, for small t one has M ^ — 2eoa^ which implies phase separation. 
However, the hole mass M decreases with increasing t and can even become nega- 



tive 



41 



Hence, for sufficiently large t, one indeed obtains a doped antiferromagnet 
which is stable against phase separation. It should be noted that a negative rest 
mass M is not at all problematical from a theoretical point of view. It only means 
that the doped antiferromagnet has a lower energy than the undoped system. Still, 
since particle number is conserved, the system cannot lower its energy by creating 
fermions. 



4.2 Fermionic Contribution to the Energy 

In this subsection we compute the fermionic contribution to the energy, keeping 
the parameters and cf of the spiral fixed. For the moment, we ignore the 4- 
fermion couplings. The considerations of this paper are valid only if the 4-fermion 
couplings are weak and can be treated in perturbation theory. Furthermore, we may 
neglect the hole-two-magnon vertices proportional to A^i and N2 which involve two 
spatial derivatives and are thus of higher order than the hole-one-magnon vertex 
proportional to A. The Euclidean action of eq. ()2.16p then gives rise to the fermion 
Hamiltonian 



+ a/^(Divp{t^^^/ + D2^{^D,^{) + Ai<ijfjvl¥_, + a/vl>/t^|v^/ j] , (4.4) 

with the covariant derivative 

A*±(a;) = di^^iix) ± ivf{x)<ili{x). (4.5) 

Here \E'{^(a;) and "^{{x) are creation and annihilation operators (not Grassmann 
numbers) for fermions of flavor f = a, f3 and spin s = +,— (parallel or antiparallel to 
the local staggered magnetization), which obey canonical anticommutation relations. 
As before, aa = 1 and ct/j = — 1. The above Hamiltonian is invariant against time- 
independent U{l)s gauge transformations 

^£(x)' = exp(±za(a;))^£(x), 

vfix)' = vf{x) - dia{x), vfix)' = vf{x) exp(±2m(x)). (4.6) 

Here we consider holes propagating in the background of a conflguration with a 
spiral in the staggered magnetization. Based on the considerations in the appendix, 
we can then limit ourselves to 

vfixy = c^, vtixy = c,eR. (4.7) 
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Hence, after the gauge transformation, the fermions experience a constant composite 
vector field Vi{x)'. The Hamiltonian can then be diagonahzed by going to momentum 
space. Since magnon exchange does not mix the flavors, the Hamiltonian can be 
considered separately for / = a and f = P, but it still mixes spin s = + with s = —. 
The single-particle Hamiltonian for holes with spatial momentum p = {pi,p2) takes 
the form 



Hf{p) 



fe+c?)^ , ^ (pi+cf)(p2+cg) 



(4.8) 

The hole-one-magnon vertex proportional to A mixes the spin s = + and s = — 
states and provides a potential mechanism to stabilize a spiral phase. The diago- 
nalization of the above Hamiltonian yields 

Ei(p) = M + ^l±^ + a,^-^El±M 
2M' ^ M" 




3 I 3\ 2 

M" 



A2(ci + a/C2)2. (4.9) 



In particular, mixing via the A vertex lowers the energy eL and raises the energy 
e!^. It should be noted that in this case the index ± no longer refers to the spin 
orientation. Indeed, the eigenvectors corresponding to are linear combinations 
of both spins. The minimum of the energy is located at p = for which 

^i(O) = M + f^ + af^^ ±A\c^ + (4.10) 

Since c? does not affect the magnon contribution to the energy density, we fix it 
by minimizing eL{0) which implies cf = = 0. According to eq. ()3.7j) this implies 
that ^0 = f , i-e. the spiral is along a great circle on the sphere S"^. By repeating 
the whole calculation including terms of (9((cf)^), we have verified a posteriori that 
putting cf = is indeed justified. At present we cannot exclude that there might 
be a minimum with a lower energy for large values of the cf. Investigating this 
issue would require a somewhat tedious numerical calculation which we have not 
yet performed. For cf = = the energies of eq. ()4.9|l reduce to 

ELiP) =M + ^ + cTf^ ± A|ci + afC2\. (4.11) 
Consequently, the filled hole pockets P£ (with M' < M") are ellipses determined by 

^ + a/-^ = Tl (4.12) 
2M' ^ M" ^' ^ ' 

where is the kinetic energy of a hole in the pocket at the Fermi surface. The 
area of an occupied hole pocket determines the fermion density as 

1 /" 1 M'M" 

(27r)2 Jp£ 27r ^/M"^ - M'^ 
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The kinetic energy density of a filled pocket is given by 



The total density of fermions of all flavors is 



n = n 



" + + + n^! = ^Mes (T° + T^ + T^ + T^), (4. 15) 

/TT 

and the total energy density of the holes is 

e/. = e: + e^+ej + (4.16) 

with 

ei = iM±A\ci + afC2\)ni + ti. (4.17) 

The filling of the various hole pockets is controlled by the parameters t£ which must 
be varied in order to minimize the energy while keeping the total density of holes 
fixed. We thus introduce 

S = eh-lin, (4.18) 
where /i is a Lagrange multiplier that fixes the density, and we demand 
dS 1 

dT^ = ± ^1^1 + ^/C2| +Tl-fi) = 0. (4.19) 

One may wonder if the density of holes of flavor a and j3 should not be fixed sepa- 
rately. After all, flavor is conserved due to the accidental U{1)f symmetry. While 
the U{l)p symmetry arises for the leading terms in the effective action, it is not 
present at the microscopic level. Although they enter the effective theory only at 
higher orders, there are physical processes that can change flavor. Despite the fact 
that such processes are rare, it would hence be inappropriate to fix the fermion 
numbers of different flavors separately. 



4.3 Four Populated Hole Pockets: Homogeneous Phase 

We will now populate the various hole pockets with fermions. First, we keep the 
configuration of the staggered magnetization fixed and we vary the in order to 
minimize the energy of the fermions. Then we also vary the parameters q of the 
magnon field in order to minimize the total energy. One must distinguish various 
cases, depending on how many hole pockets are populated with fermions. In this 
subsection, we consider the case of populating all four hole pockets (i.e. with both 
flavors f = a,l3 and with both energy indices ±). In this case, eq. ()4.19p implies 

vrn vm . , , , , 
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Figure 2: The homogeneous phase with constant staggered magnetization. 



The total energy density then takes the form 

e = eo + em + eh = eo + 2p,{cl + cl) + el + e'l+e1 + e^_ 

TITX^ 1 

= to + 2ps{c\ + cl) + Mn + — M,^K\cl + cl). (4.21) 

4Mefr TT 

For 2Txps > McffA^ the energy is minimized for Q = and the configuration is thus 
homogeneous. The total energy density in the four-pocket case is then given by 

" = ^» + *'"+4Mi' 

The homogeneous phase is shown in figure 2. For 27rps < MeffA^, on the other 
hand, the energy is not bounded from below. In this case, cf + seems to grow 
without bound. However, according to eq. ()4.20p this would lead to Tjf < which 
is physically meaningless. What really happens is that one pocket gets completely 
emptied and we must thus turn to the three-pocket case. 



4.4 Three Populated Hole Pockets: 45 Degrees Spiral 

We now populate only three pockets with holes: the two pockets with the lower 
energies i?" and as well as the pocket with the higher energy E'^. Of course, 
alternatively one could also fill the /?+-pocket. We now obtain 

n = n^ + n^+n(! = ^Meff(T^ + r^ + rf), = + el + e^_, (4.23) 
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such that eq. ()4.19|) yields 



2?™ A, , 
u = M H — Ci — Co , 

™ 27m , , , A , , 

T, = — TT A Ci + Co Ci — Co , 

+ 3^^^ II 21 3I 1 21, 

^rv 27rn , , , A , , 

T = — TT h A Ci + Co Ci — Co , 

a 27m 2A, , , , 

The total energy density then takes the form 

e = eo + e^ + eh = eo + 2ps{cl + cl) + el + e'^ + e'^ 

/ On / A, A 7rn^ 

= eo + 2p,(cf + c^)+ M--|ci-C2|U + 



V 3' ^ ^7 3Meff 
- ;^MeffA2(c? + C1C2 + cl). (4.25) 

For 2nps > M^ffA^ the energy density is bounded from below and its minimum is 
located at Ci = — C2 with 

II II 

This represents a spiral in the staggered magnetization oriented along a lattice di- 
agonal — a 45 degrees spiral. When one occupies the instead of the a+-pocket, 
one finds Ci = C2, i.e. the spiral is then oriented in the orthogonal diagonal direction. 
According to eq. ()A.19jl in the appendix, the spiral pitch (i.e. the wave number of 
the spiral) is given by 

c3c3 + c.c.^2,/c?+cH^^^^^^. (4.27) 
The resulting energy density in the three-pocket case takes the form 

TV ( ^ 1 MeffA^ ^ 2 



eg = Co + Mn + 1 n\ (4.28) 

The 45 degrees spiral is illustrated in figure 3. It is energetically less favorable than 
the homogeneous phase because €3 > for 271 ps > M^sA'^. Only for 2nps = MefjA^ 
both phases cost the same energy, i.e. €3 = €4. For 27ips < MefjA^ the energy density 
of eg. ()4.25|1 is unbounded from below and cl + cl again seems to grow without bound. 
This, however, would lead to T° < which is unphysical. In fact, the a+-pocket is 
completely emptied and we are thus led to investigate the two-pocket case. 
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Figure 3: A 45 degrees spiral oriented along a lattice diagonal. 



4.5 Two Populated Hole Pockets: Zero Degree Spiral 

We now populate only two pockets with holes. These are necessarily the pockets 
with the lower energies E'^ and E^. In this case we have 

n = n° + = ^M,s{T^ + tH), = e'l + e^_, (4.29) 
Ztt 

and thus eq. ()4.19p now implies 

Tin A , , , , 

l^ = ^+J^- 9 (1^1 + + |Cl - C2I), 

T^n A , , 

T° = -— + -(|Ci + C2|-|Ci-C2|), 

T^ = ^ + ^{\ci-C2\- |ci + C2I). (4.30) 

The total energy density then takes the form 

e = 60 + em + eh = €0 + 2ps{cl + cl) + e'l + 

= eo + 2ps{cl + 4)+ [M- -(|ci + C2I + |ci - C2I) n + 



1 A^ 

- — Meff— (|Ci+C2|-|Ci-C2|)l (4.31) 
ZTT 4 

For 27rp, > ^M^sA^ the energy density is bounded from below and has its minimum 
at 

, , A , , , , , , A , , 

ci = - — n, C2 = 0, or ci = 0, C2 = - — n, (4.32) 
4ps 4ps 
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Figure 4: A zero degree spiral oriented along a lattice axis. 



which corresponds to a spiral in the staggered magnetization oriented along a lattice 
axis — a zero degree spiral. Again, according to eq. ()A.19|) . the wave number of the 
spiral is given by 



The zero degree spiral is more stable than the homogeneous phase if €2 < 64, which 
is the case for 27r ps < MefjA^. This also clarifies the instability of the homogeneous 
phase in this parameter regime: it simply turns into the zero degree spiral. We have 
verified explicitly that the assumption cf = is justified by repeating the whole 
calculation including terms up to C((cf)^). Indeed, cf = turns out to be a stable 
minimum for 2ttPs < MegA^. The zero degree spiral is illustrated in figure 4. 

4.6 One Populated Hole Pocket: 45 Degrees Spiral 

Finally, let us populate only one hole pocket, say the states with energy E^. Of 
course, alternatively one could also occupy the /3_-pocket. One now obtains 




(4.33) 



and the resulting energy density in the two-pocket case takes the form 




(4.34) 



T 



27m 



(4.35) 



Meff 
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The total energy density then takes the form 



eo + Cm + e/i = eo + 2ps{cl + cl) + e°_ 
eo + 2ps{cl + cl) + (M - A|ci + CaDn + 



Tin 



eff 



which is minimized for Ci = C2 with 



|Cl| 



A 

4p. 



-n. 



(4.36) 



(4.37) 



This again represents a 45 degrees spiral in the staggered magnetization oriented 
along a lattice diagonal. When one occupies the /5_-pocket one finds ci = — C2, i.e., 
as in the three-pocket case, the a- and /3-spirals are oriented in orthogonal diagonal 
directions. Note that the three-pocket 45 degrees spiral with an occupied a+-pocket 
is oriented in a direction orthogonal to the one of the one-pocket 45 degrees spiral 
with an occupied a_-pocket. In the one-pocket case, the spiral pitch is given by 



k 



cjc^ + CiCi 



cl + cl 



An 



(4.38) 



and the corresponding energy density takes the form 



ei = eo + Mn + — — 



n 



(4.39) 



The 45 degrees spiral is energetically more favorable than the zero degree spiral 
if ei < €2, which is the case for 27rps < |MefrA^. Again, by repeating the whole 
calculation including terms up to (9((c^)^), we have verified explicitly that for 27ips < 
MefjA^ the assumption of cf = is justified because it corresponds to a stable 
minimum. 



However, as already pointed out in 24, 27, 3^ i^, 0]' negative coefficient of 
the term proportional to n"^ (i.e. a negative compressibility) leads to an instability of 
the 45 degrees spiral against increasing the local fermion density by decreasing the 
wavelength 27i/k of the spiral. To see this, let us consider a one-pocket 45 degrees 
spiral filling a sub-volume V and leaving an undoped antiferromagnet behind in 
the remaining volume V — V. The energy density of this configuration with an 
inhomogeneous fermion density is given by 



V 



V -V 



Mn + 



7r 



A^ 



eo + 



M 



TT 



TT 



A^\ /iV 
off '^sj \V' 

2\ 



Meff Aps 
21 y, 



n 



YL 

V 



A' 



M, 



eflf 



ApsJ V'\ V 



N 



(4.40) 
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Since the coefficient of the term proportional to 1/V is negative for 2t{ ps < |MeffA^, 
the system can minimize its energy by shrinking V and thus increasing the fermion 
density n in the region of the 45 degrees spiral. The one-pocket 45 degrees spiral is 
thus unstable against the formation of inhomogeneities in the fermion density. Our 
basic assumption that the system is homogeneous is then not satisfied. 

4.7 Symmetry Breaking Pattern in the Spiral Phase 

As we have seen, for intermediate values of Ps a spiral phase replaces the phase 
with homogeneous staggered magnetization. In the homogeneous phase the SU{2)s 
spin symmetry is spontaneously broken to its U{l)s subgroup. Due to antiferromag- 
netism, the displacement by one lattice spacing is also spontaneously broken. In a 
spiral phase, on the other hand, a U{l)s spin rotation that leaves the configuration 
invariant no longer exists. Hence, U{l)s is then also spontaneously broken and one 
may expect an additional massless Goldstone boson. Furthermore, due to the finite 
wavelength 2TT/k of the spiral, translation invariance (not only by one lattice spac- 
ing) is also spontaneously broken. Only the translations by an integer multiple of the 
spiral wavelength remain unbroken. Due to the spontaneously broken translation 
symmetry one expects yet another massless Goldstone boson — a spiral phonon (or 
helimagnon) — which corresponds to vibrations of the spiral. In order to correctly 
count the number of Goldstone bosons one must notice that while both U{l)s and 
translation invariance are individually spontaneously broken, there is a combination 
of these two symmetries that remains unbroken. In particular, any translation of 
the spiral configuration can be compensated by an appropriate C/(l)s spin rotation. 
Consequently, there are not two but there is only one additional Goldstone boson 
— the spiral phonon. Besides this additional massless boson, in the spiral config- 
uration there are still two magnons — a transverse and a longitudinal one. The 
transverse magnon corresponds to spin fiuctuations out of the spiral plane, while 
the longitudinal magnon represents in-plane fiuctuations. 

It should be noted that the wavelength 27T/k oc 1/n of the spiral represents a 
new length scale in the problem. In particular, for very small fermion density n 
this length scale is arbitrarily long. In that case, the spiral phonon will have very 
little effect on the dynamics of the holes. When the fermion density increases, the 
spiral wavelength shrinks and the spiral phonon becomes more important. In partic- 
ular, besides magnon exchange the exchange of spiral phonons may then contribute 
significantly to the long-range interactions between the holes. 

Since the symmetry is no longer broken just from SU{2)s to U{l)s, one may 
wonder whether the vector e{x) living in the coset space SU{2)s/U{l)s = S"^ is still 
an appropriate low-energy degree of freedom in the spiral phase. Fortunately, this 
is indeed the case. Since the spiral phonon arises as an additional Goldstone boson, 
one may also wonder whether a new field must be added to the effective Lagrangian. 
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This is not necessary because, just like the two magnons, the spiral phonon arises 
as a fluctuation of the staggered magnetization vector e{x) (in this case around the 
spiral configuration). Only if one would construct another effective theory valid only 
at extremely long length scales much larger than the wavelength of the spiral, the 
basic fields of the theory would have to be redefined. In such an effective theory 
the spiral itself would be a short distance phenomenon and the spiral phonon would 
appear explicitly as an independent degree of freedom. The corresponding physics 
at extremely low energies is still contained in our effective theory which is valid at 
length scales both longer and shorter than the wavelength of the spiral. 



5 Inclusion of 4-Fermion Couplings in 
Perturbation Theory 

In this section the 4-fermion contact interactions are incorporated in perturbation 
theory. Depending on the microscopic system in question, the 4-fermion couplings 
may or may not be small. If they are large, the result of the perturbative calculation 
should not be trusted. In that case, one could still perform a variational calculation. 
In this work we limit ourselves to first order perturbation theory. We will distinguish 
four cases: the homogeneous phase, the three-pocket 45 degrees spiral, the zero 
degree spiral, and the one-pocket 45 degrees spiral. Finally, depending on the values 
of the low-energy constants, we determine which phase is energetically favorable. 



5.1 Four-Pocket Case: Homogeneous Phase 

Let us first consider the homogeneous phase. The perturbation of the Hamiltonian 
due to the leading 4-fermion contact terms of eq. ()2.18|) is given by 

AH = J £x J2 [^(^f^:*-l^^. + *f^f^^t^^J 

+ Gs^f + Gg^f J . (5.1) 



s=+, 



It should be noted that ^{^(x) and "^{{x) again are fermion creation and annihilation 
operators (and not Grassmann numbers). The term proportional to G4 is of higher 
order and hence need not be considered here. In the homogeneous phase the fermion 
density is equally distributed among the two spin orientations and the two flavors 
such that 

(^^t^a^ = (^°t^^^ _ (^f = _ ^_ (5 2) 

The brackets denote expectation values in the unperturbed state determined in 
section 4.3. Since the fermions are uncorrelated, for f ^ f 01 s ^ s' one has 

(^/tv^/vl/Z'tvI;/;) = (^/t^/)(^/;t^/y (5.3) 
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Taking the 4-fermion contact terms into account in first order perturbation theory, 
the total energy density of eq. ()4.22|) receives an additional contribution and now 
reads 

e, = eo + Mn + -— + -(d + + G^)n\ (5.4) 

4Mefr O 



5.2 Three-Pocket Case: 45 Degrees Spiral 

For = the eigenvectors of the single-particle Hamiltonian of eq. ()4.8|1 correspond- 
ing to the energy eigenvalues Ej^{p) are given by 

= ^{^L±^{) ^1 = ^{^{^^i). (5.5) 
V 2 v 2 

Inserting this expression in eq. ()5.1|) allows us to evaluate the expectation value ( AiJ) 
in the unperturbed states determined in section 4. In the three-pocket case the states 
with energies E'^{p), E^{p), as well as E'^{p) (or alternatively E^{p)), and with p 
inside the respective hole pocket are occupied and 



As a result the energy density of eq. ()4.28p turns into 
eg = eo + Mn + 1 — 

+ I [{Gi + G, + G,)4np., - GiMef^A^] ^vrp. - MefrAj 
4 (67rps - McffA^)^ 



5.3 Two-Pocket Case: Zero Degree Spiral 

In the zero degree spiral only the states with energy E"{p) and E^{'p) with p inside 
the respective hole pocket p!_ are occupied and hence 

($-t5«^ = ($/^t^/3^ ^ ($?$°) = = 0. (5.8) 

As a result the energy density of eq. ()4.34|) turns into 
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5.4 One-Pocket Case: 45 Degrees Spiral 



In the one-pocket case only the states with energy E'^{p) (or alternatively with 
E^{p)) and with p inside the corresponding hole pocket are occupied such that 

(^"t^a^ = ($?$°) = = {^>^W_) = 0. (5.10) 

In this case the 4-fermion terms do not contribute to the energy density which thus 
maintains the form of eq. ()4.39|) . i.e. 



5.5 Stability Ranges of Various Phases 

Let us summarize the results of the previous subsections. The energy densities of 
the various phases take the form 

ei = to + Mn + ]^Kin^. (5.12) 

According to eqs. ()5.1H) . (j5.9|l . (j5.7|l . and ()5.4|1 . the compressibilities are given by 
27r A2 



271 ( 1 MeffA2 

'^S = 7771— 1-77: 



3Meff V 2 67rp, - MeffA2 
+ \ [{G, + G, + G3)47rp. - GiMeffA^] ^^^^ ~ 



= ^ + i(^i + ^2 + G3). (5.13) 

The compressibilities Hi as functions of M^gA"^ /2'n'Ps are shown in figure 5. For large 
values of ps, spiral phases cost a large amount of magnetic energy and a homogeneous 
phase is more stable. The energy density of the homogeneous phase is lower than 
the one of the spiral phases for K4 < /ti, K2, K3, which is the case for 

, , ,2 ^ M^fT-A^ r Four-pocket case: ,x 

MeffA^ + Gi— £S — <2'Kps < TT 1 (5.14) 

27r 1^ Homogeneous phase. 

At 0{Gi) the three-pocket 45 degrees spiral always costs more energy than the other 
phases, except at the specific point 

71 ^ A 2 ^ M^ffA'^ f Three-pocket case: -.^x 
2ttPs = M,sA^ + Gi^^ 1 -1 5.15 

27r I 45 degrees spiral. 
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Homogeneous phase 



Zero degree spiral 



2 M,jfA^/27ip, 
Inhomogeneous phase 



Figure 5: The compressibilities Ki (in the presence of weak repulsive 4-fermion cou- 
plings) as functions of M gffA^ /27r Ps determine the stability ranges of the various 
phases. A homogeneous phase, a zero degree spiral, or an inhomogeneous phase are 
energetically favorable, for large, intermediate, and small values of ps, respectively. 
Only at the two isolated points separating the three regimes a 45 degrees spiral, either 
with one or with three populated hole pockets, is degenerate with the other phases. 
The three- and four-pocket cases become unstable at Mf^gt^ jl'Kps — 1, and the one- 
and two-pocket cases become unstable at M^gb3- j^'Kps — 2. 



For this particular value of ps the three-pocket 45 degrees spiral costs the same 
energy as the homogeneous phase and the zero degree spiral. For larger values 
of the 4-fermion couplings (not accessible to first order perturbation theory) it is 
conceivable that the 45 degrees spiral becomes energetically favorable. This could be 
investigated, for example, using a variational calculation. For intermediate values of 
Ps the homogeneous phase becomes unstable against the formation of a zero degree 
spiral, i.e. k.2 < ka, and we find 

Ve.A^ + (G2 + Gs)^< 2nps < MeffA^ + G,^^^^" ^ Two-pocket case: 



2 ^ csy _ - en 27r 1 Zero degree spiral. 

(5.16) 

For smaller values of ps the one-pocket 45 degrees spiral becomes energetically more 
favorable than the zero degree spiral, i.e. Ki < K2, but the one-pocket 45 degrees 
spiral exists only at the isolated point 

= + (G, + G3)^ I (5.17) 

2 Stt I 45 degrees spiral. ^ ^ 
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It may seem that a repulsive 4-fermion interaction G2 + > stabilizes the 45 
degrees spiral, at least in the narrow parameter range down to 2nps = ^M^sA^. 
However, although we have not yet identified the nature of the inhomogeneous phase, 
we expect that it will be energetically more favorable, such that 

2.P. < + (G, + G,)Mk^ I ^rr^ (6.18) 



2 ^ Stt \ yet unidentified nature. 

At least at 0{Gi), the 45 degrees spiral cannot be realized in this regime, because 
for Gj = it only exists at the isolated point 27rps = ^MeffA^ and first order 
perturbation theory uses just the unperturbed wave function. Definitely, the one- 
pocket 45 degrees spiral becomes unstable when ki < 0, i.e. for 

r. ^ 1 71 A 2 f Instability of the 45 degrees spiral 
27rps < -MeffA < ■ ^ r X- r • r -x- (5-19) 

2 [ agamst formation of mhomogeneities. ^ ^ 

It should be pointed out again that these results apply only if the 4-fcrmion contact 
interactions are weak. If the 4-fermion couplings are strong, the approach of filling 
pockets with weakly interacting holes is not applicable and the true ground state of 
the system may be different. Even if the 4-fermion couplings are small, the present 
result does not necessarily reveal the complete nature of the true ground state. In 
particular, the configurations of the staggered magnetization were restricted to those 
that imply a homogeneous composite vector field Vi{xy, i.e. to spirals or to config- 
urations that are homogeneous themselves. For example, the double spiral, which 
implies an inhomogeneous Vi{xy, was not taken into account. Hence, we cannot 
exclude that the phases that we found may still be unstable against developing in- 
homogeneities, at least in certain parameter regions. This shall be explored in the 
future. 



6 Reduction of the Staggered Magnetization 
upon Doping 

Until now we have considered the Lagrangian or the Hamiltonian of the effective 
theory. Now we will see that other quantities can be constructed in a similar way. 
An observable of particular interest is the staggered magnetization which is the 
order parameter for the spontaneous breakdown of the SU{2)s symmetry. As the 
antiferromagnet is doped, the staggered magnetization is reduced until the SU{2)s 
symmetry is restored. In the actual materials this happens at relatively small doping, 
before high-temperature superconductivity sets in. 

In the undoped antiferromagnet the local staggered magnetization is given by 
Ms{x) = M.sS{x). The low-energy parameter M.s is the length of the staggered 
magnetization vector. For the Heisenberg model (or equivalently for the t-J model 
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at half-filling) this parameter has been determined with high precision in quantum 
Monte Carlo simulations resulting in Aig = 0.3074(2)/a^ jT^, E^- In a doped 
antiferromagnet the staggered magnetization receives additional contributions from 
the fermions such that 

M,{x) = Ms-mY^ iljp{x)ip{{x) e{x). (6.1) 

s = +,- 

Here m is another low-energy parameter which determines the reduction of the 
staggered magnetization upon doping. It should be noted that there are further 
contributions to Ms{x) which include derivatives or contain more than two fermion 
fields. All these terms are of higher order and will be neglected here. Since both 
the homogeneous and the spiral phases have a constant fermion density we can use 

Y{^M)=n, (6.2) 



such that 

Ms{n) = Ms-mn. (6.3) 

The higher-order terms that we have neglected correct this equation at Oln"^). A 
rough estimate of the critical density at which the SU (2)^ symmetry gets restored is 
ric ~ M-s/m. It would be interesting to determine the value of m for the t-J model, 
which may be feasible in quantum Monte Carlo calculations. 



7 Conclusions and Outlook 



In this paper we have used a systematic effective field theory for antiferromagnetic 
magnons and holes in order to investigate the propagation of holes in the background 
of a staggered magnetization field. We have limited ourselves to configurations of 
the staggered magnetization that are either homogeneous themselves or that gen- 
erate a constant composite vector field Vi{x)' for the fermions. In both cases, the 
resulting fermion density is homogeneous. Our calculations also rely on the as- 
sumption that remnant 4-fermion contact interactions between the holes are weak 
and can be treated in perturbation theory. This may or may not be the case for a 
concrete magnetic material. We like to emphasize again that the effective theory 
is universal and makes model-independent predictions for a large class of magnetic 
systems. Material-specific properties enter the effective theory through the values of 
low-energy parameters such as the spin stiffness ps- For large values of ps distortions 
in the staggered magnetization cost a large amount of energy and a homogeneous 
phase is energetically favored. In that case, all four hole pockets are equally pop- 
ulated with doped holes. For smaller values of p^, on the other hand, the doped 
holes can gain energy from a spiral in the staggered magnetization. For intermedi- 
ate values of ps a zero degree spiral, in which only two hole pockets are populated. 
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is most stable, while a 45 degrees spiral (with either one or three populated hole 
pockets) can exist only at two specific isolated values of ps- It is conceivable that the 
45 degrees spirals may be stabilized for larger values of the 4-fermion couplings Gi 
(inaccessible to first order perturbation theory). For small values of ps a yet uniden- 
tified inhomogeneous phase is energetically favored. The reduction of the staggered 
magnetization upon doping is again controlled by a low-energy parameter m whose 
value depends on the material in question. 

The results of our investigation provide a basis for further studies that naturally 
suggest themselves. First, it would be interesting to investigate if there is a stable 
minimum of the energy for large values of cf. This requires a straightforward but 
somewhat tedious numerical calculation which we have not yet performed. Next, 
instead of using first order perturbation theory, one may want to include the 4- 
fermion couplings in a variational calculation. This would eliminate the assumption 
that these couplings are small. One should also analyze the stability of the homoge- 
neous and spiral phases against developing inhomogeneities in the fermion density. 
For example, it is interesting to ask if a double spiral is energetically more favorable 
than the spiral phases considered here. If this were the case, it is conceivable that a 
lightly doped antiferromagnet without impurities is an insulator. The real materials 
are indeed insulators, but they contain impurities on which the doped holes may get 
localized. 

On the other hand, if the homogeneous phase or the simple spirals considered here 
turn out to be more stable than, for example, a double spiral, the effects of magnon 
exchange between doped holes would be interesting to study in detail. In particular, 
using the effective theory, we have shown that the one-magnon exchange potential 
between two isolated holes gives rise to binding with (i-wave characteristics 71, 7^. 



In a spiral phase the exchange of spiral phonons (i.e. helimagnons) may also be an 
important dynamical mechanism for the binding of hole pairs. Depending on the size 
of the pairs and on their density, hole pair formation may lead to Bose-Einstein con- 
densation or to magnon-mediated BCS-type superconductivity coexisting with an- 
tiferromagnetism. While coexistence of superconductivity and antiferromagnetism 
is not observed in the high-Tc cuprates, this may be due to impurities created by 
doping on which holes may get localized, thus preventing superconductivity. A clean 
system like the t-J model, on the other hand, may superconduct already at small 
doping within the antiferromagnetic phase, although T^. may then be rather small. 
The effective theory provides us with a tool that allows us to address such questions. 
If there is superconductivity already within the antiferromagnetic phase in a clean 
system, the corresponding mechanism responsible for superconductivity may persist 
in the high-Tc cuprates, despite the fact that in the real materials at small doping 
superconductivity may be prevented by the localization of holes on impurities. Of 
course, once antiferromagnetism is destroyed, magnons and spiral phonons no longer 
exist as massless excitations. However, antiferromagnetic correlations, although only 
of finite range, still exist in high-Tc superconductors and may still play an important 
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role as relevant low-energy degrees of freedom. This is similar to nuclear physics 
where the pion is not exactly massless (in that case due to explicit chiral symmetry 
breaking) but is certainly relevant at low energies. Still, whether magnon-mediated 
binding between holes may be responsible for high-temperature superconductivity 
remains a difficult question that may or may not be within the applicability range of 
the effective theory. For lightly doped cuprates the low-energy effective field theory 
is applicable. It yields reliable and interesting results and should be pursued further. 
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A Most General Configuration with a 
Homogeneous Composite Vector Field 

In this appendix we show that the most general configuration with constant Vi{x)' 
corresponds to a spiral in the staggered magnetization. According to eg. 1)3.11) we 
need to consider 9{x), (p{x), and a{x) such that 

9(x) 

'^t^'A' — "^^^(3;) exp(±2ia(x)) 

= ^(sin6'(a;)9jV9(a;) ± Z(9j6'(a;)) exp(=Fi((y9(2;) — 2a;(x))) = , (A.l) 

with and cf being constant. Introducing x(a^) = V^(a^) — 2a (x) one obtains 
cos 9{x)dia{x) = sin^ ^ dix{x) — cf, 

tan6{x){dixix) — 2cf) ±idi6{x) = 2c^ exp{±ix{x)) = 2cj exp(±i(x(x) +uji)). 

(A.2) 

Here we have put 

cf = Ciexp{±iuJi), with q G M, ^^i ^ [-^,^]- (A. 3) 
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By an appropriate constant gauge transformation, one can always put uJi = 0. From 
eq. ()A.2|) one infers 



(A.4) 



tane{x){dixix) -2c-) =2qcos(x(x) +uji), 
di6{x) = 2ci sin(x(a;) + tOi). 

Demanding did20{x) = d2di6{x) leads to the constraints 

= a, Co"! = Co'2 = 0, 



(A.5) 



which imply 



diXi^) = 2cf [1 + a cot 6'(x) cosx(x)], 
di9{x) = 2c^asmx{x). 



(A.6) 



These equations can be satisfied only if x and 6 are functions of 2; = cfxj, i.e. if they 
are plane waves with 



dzXi^) = 2[1 + a cot 6{z) cos x{^)]y 
dz6{z) = 2asinx(2;). 



(A.7) 



Since both x and 9 depend on z only, we may also consider % as a function of 9 
which then leads to 



dexiO) 



This can be cast into the form 



a sin x{9) 



1 



+ cot6'cotx(6'). 



-dgcosx{9) = — h cot ^ cos x(^), 
a 



and is solved by 



cosx(6') 



cos 9 — X 
a sin 9 ' 



(A.8) 



(A.9) 



(A.IO) 



where A is an integration constant. Inserting this result in eq. ()A.7|) one obtains 



d,cos9{z) = ±2v/a2(l - cos"^ 9{z)) - {cos 9{z) - A)2, 
2{cos 9{z) - A) 



sin^ 9{z) 



(A.ll) 



The equation for dz^p{z) results by combining eq. ()A.2|) with eq. ()A.7|) and eq. ()A.10|) . 
The above equation for cos^^(2;) can be integrated in a straightforward manner and 
one obtains 



cos6'(2;) 



l + a2 
1 



A + aVl + a2 - A2 cos 2Vl + a?(z - 



cos 1] + asmi] cos ( 2\^l + a^{z — Zq] 
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(A.12) 



Here we have expressed the integration constant as 

A = Vl + a^cosr/. (A.13) 

— * 

It will turn out that r] is the angle between the direction j perpendicular to the 
spiral plane and the 3-axis. Furthermore, one obtains 

sin 9 (z) cos {ip{z) — Lp{zo)) = = sin r/ — a cos t] cos (2\/ 1 + a^^z — zq)) , 

sin e{z) sm{<^{z) - <^{zo)) = " sin (2VI + a^{z - Zo)) , (A.14) 
and thus 

/X / a sin (2VI + 0^(2; — 2:0)) \ /. X 

^(z) - ^(^0) = atan ^ ' . (A.15) 

\ sm ?7 — a cos ?7 cos (^iv 1 + 0, (-2 — ^o) j / 

Differentiating this equation with respect to z, it is straightforward to show that 
eq. ()A.ll|) is indeed satisfied. 

We now form the scalar product of the unit- vector 

e{z) = (sin6'(2) cosv9(2;), sin6'(z) siny9(2;), cos^^(2;)), (A. 16) 

describing the staggered magnetization, with the unit-vector 

j = (sin?7 cos v5(-Zo), sin?7 sin (y9(zo), COS?]), (A. 17) 

which yields 

e{z) ■ j = sin6'(z) sin r] {cos {p{z) cos ip{zQ) + sirnp^z) sin ip^zo)) + cos6'(z) cos 77 

= sin 6{z) sin r] cos(v?(z) — fjzo)) + cos 9{z) cos r] = —= , (A. 18) 

V 1 + 

i.e. a constant (2;-independent) scalar product. This finally proves that e{z) indeed 
describes a spiral in a plane perpendicular to j. Replacing z = c^Xi and using 
eq. ()A.5jl . the wave number of the spiral (i.e. the spiral pitch) can be identified as 



k = 2Vl + a^^cfcf = 2^44 + dd. (A. 19) 

The particular spiral configuration considered in eq. ()3.4|) has 77 = such that 

cos6'(z) = —^=^= = cos 6*0, 
Vl + 

sin 6{z) sin{(f{z) — fizo)) = —^=^= sin [2^/1 + a?{z — zq) 

yl + a? ^ 

if{z) - if^zo) = -2VI + = kiXi. (A.20) 

This is indeed consistent with eq. ()3.8j) which implies 

a = -tan^o, h = hr- (A.21) 

cost^o 
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